RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D, No. 9, September 1965 



Magneto-Ionic Propagation in Inhomogeneous 
Media. Part I Transverse Propagation 

B. K. Banerjea 
Institute of Armament Technology, Poona, India 

(Received March 29, 1965) 

Propagation of electromagnetic waves in an inhomogeneous magneto-ionic medium has been dis- 
cussed with special reference to ionospheric propagation. This involves the study of two coupled 
linear second-order differential equations with coefficients having one essential and more than one 
regular singularity. The problem has been tackled in two stages, divided into parts I and II of this 
paper. 

In part I, dealing with the transverse propagation in a parabolic electron layer, the wave equation 
belonging to a class higher than the confluent hypergeometric or Mathieu equations has been investi- 
gated. In part II, dealing with the oblique propagation, the effect of the coupling function for the above 
class of equations has been studied, and relevant asymptotic solutions for the coupled equations 
derived. 

Mathematical expressions for the physical observables in ionospheric investigations have been 
derived from the solutions of the wave equations inside a parabolic electron layer. 

1. Introduction 

A partially ionized medium traversed by a steady magnetic field is termed a magneto-ionic 
medium; the ionosphere is such a medium, in which the electron density is a function of the vertical 
height. It can be shown from the Maxwell-Lorentz equations that the propagation of electromag- 
netic (em) waves inside such a medium is controlled by the equation 
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where mks units and symbols (see appendix) have been used. 

The standard pulse technique method used for the study of ionospheric propagation furnishes 
information about the following observables: 

(a) The interval r taken by the pulse to travel from the ground to the ion layer and back. 

(b) The fraction R 2 of energy (reflection coefficient) reflected back by the ion layer. 

(c) The state of polarization of the downcoming wave as determined by the ratio R of the polari- 
zation ellipse and the tilt angle *P of the major axis with respect to the local magnetic meridian. 

Any method developed for the study of magneto-ionic propagation of em waves must succeed 
in furnishing mathematical expressions in terms of the known parameters for these observables 
from the solutions of (1). 

The ray optics of Appleton [1927] and Hartree [1929] form a special solution of (1) with the 
assumptions: 

(i) — =2m/, where /is the frequency of the incident wave, and 
at 

(ii) N=No, a constant. 
The observed values of r are utilized in the ray optics for ascertaining the electron density N as 
a function of the equivalent height P' = \ct, a distance which light would have covered in an in- 
terval of r/2. For frequencies less than the plasma frequency of the electron layer, the reflection 
coefficient R 2 as computed by the ray optics is unity for a friction-free layer. Any divergence from 
this theoretical value is interpreted as due to frictional absorption, and appropriate values of the 
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collision frequency v are affixed accordingly. Thus two out of the four observables are utilized 
for predicting the values of N and v. The state of polarization as deduced by the ray optics is 
determined by the lowermost stratum with vanishingly small electron density. This prediction of 
the ray optics is contrary to actual experience. 
The basic difficulties in the study of (1) are: 

(a) Magneto-ionic anisotropy being of the cyclotonic class, separate wave equations describing 
the propagation of differently polarized waves cannot be set up unless N=N (ray optics), or 
w//Xk = (longitudinal propagation) or w// • k = (transverse propagation), k being the propagation 
vector. 

(b) Magneto-ionic propagation in inhomogeneous media is described by equations of a class 
higher than those studied in other branches of mathematical physics. 

The successive stages in the development of the general solution of (1) where the assumptions 
of the ray optics are not tenable can be arranged as follows: 

(i) Removal of the partial time derivatives in (1) by applying Laplace-transform with respect 
to time and expressing the nonhomogeneous terms arising therefrom in terms of the prescribed 
initial conditions. Though the plasma frequency \/Ne 2 lmeo of the electron layer and the gyro- 
frequency \e/JLoRelm\ of the electrons around the terrestrial magnetic field H P are of the same 
order as the frequencies of the incident waves, it can be shown that noncoherent scattering is not 
significant for the order of values of N, v, and w// inside the ionosphere. Hence the partial time 
derivatives in (1) can be replaced by the factor io) = 2i7rf without introducing any significant error. 

(ii) A two-dimensional Laplace-transform with respect to the space-variables normal to V/V 
(on the assumption that 7iVx V^ = 0) will reduce (1) into a pair of second-order ordinary linear 
differential equations, each containing two components of the electric vector (or its transform, if 

— is not replaced by ia>) normal to V TV. For vertical propagation of plane em waves this step is 

superfluous. In this latter case the equations controlling the propagation can be written as 

-|P+ [Uz)E x - iQL(z)E y ] = 0, 

(2) 
+ [K2(z)E y + iQL(z)E x ] = 0, 



cPE y 



dz 2 

where z is the vertical height and the functional dependence of Ki(z), £2(2:), and L(z) depends on the 
forms of N(z) and v(z). 

(iii) The pair of equations in (2), for given forms of N(z) and v(z), can be studied by one of the 
following methods: 

(a) Conversions into canonical form and application of the standard technique for the study of 
a system of equations [Ince, 1956]. 

(b) Application of perturbation methods to (2) or to equations derived from (2) and satisfied 
by suitable linear combinations of E x and E y . 

(c) Elimination of one of the dependent variables by raising the order of differentiation and 
studying the resulting fourth-order differential equation involving one dependent variable. 

If utilized properly and interpreted correctly, all the above methods will lead to identical 
results within a certain degree of accuracy. However, the coefficient functions K\(z), K 2 (z) and 
L(z) are so complicated that even for the simplest types of electron layers, the solutions of (2) 
cannot be expressed in closed forms in terms of functions whose properties have already been 
studied. Fortunately, expressions for the physical observables r, 7? 2 , £$, and 9 can be derived 
from the asymptotic solutions of (2). These solutions for \z\— *°° can always be expressed as 

E x ~ exp [i{cot-\-(t) lQ }]-\-R q exp [i{(ot + <f> 2q }], 

E v~^i q ex P [*{«* + *i f + *i ff }]+#«A ex P [iW + 2 , + yi (3) 
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where the suffix q refers to the qth magneto-ionic component, while /?<?, ^ 2 , </> 1(7 2g , and i// 1(/ 2(? 
are real quantities expressed in terms of the layer parameters, pulsatance ci> = 277/of the incident 
wave, and the geomagnetic characteristics of the place of observation. 

For each mode of propagation, the observables r, /? 2 , 0i 9 and ^ are related to the above asymp- 
totic solution as 

2^ 
t = (<K~<M/o>, R 2 = [R Q ]\ M=M W tan 2^^^-sin i// lg , (4) 

where ^F is the tilt angle of the polarization ellipse traced out by the magnetic vector of the reflected 
echo. 

For the transverse propagation the coupling coefficient QL(z) in (2) vanishes for any type of 
electron layer and separate equations describing the propagation of two linearly polarized waves 
can be set up. The transverse mode with the electric vector parallel to H e is generally termed as 
the ordinary or o-mode, while the other termed as extraordinary or e-mode. For the longitudinal 
modes K\{z) = K 2 (z) = K(z) for any type of electron layer, and (2) reduces to 

j^ (E x ± iE v ) + [K(z) + L{z)\E x ± iE y ) = 0. (5) 

These equations describe the propagation of two circularly polarized magneto-ionic components. 

The effect of the steady magnetic field H e on the propagation of the transverse o-mode or either 
of the longitudinal modes is nil or trivial. Hence these degenerate cases of propagation belong 
rather to the ionic than to the magneto-ionic class. Ionic propagation in linear, parabolic, and 
exponential electron layers has been studied by Epstein [1930; Hartree, 1931; Wilkes, 1940; Ryd- 
beck, 1942; Wait, 1962], and others. Solutions of the wave equations for these cases can be ex- 
pressed in terms of the functions of the confluent hypergeometric group, viz, Bessel or Weber 
functions. 

The only remaining mode of these degenerate cases is the transverse e-mode whose propaga- 
tion is characteristically influenced by the magnetic field H<>. The propagation of this mode for 
the type of electron layers mentioned above is described by equations having, besides the essential 
singularity at infinity, more than one regular singularity in the finite region. This mode represents 
the simplest type of proper magneto-ionic propagation and has not been studied by any worker as 
yet. In part I of this paper the propagation of this transverse e-mode for a parabolic electron layer 
has been investigated in detail, and theoretical expressions for the relevant observables, viz, r 
and R 2 , have been derived in terms of the layer parameters. 

As regards the general propagation with V/V 4= 0, w H X k =t= 0, oqh ' k 4= 0, the functions Ki(z), 
K 2 (z), and L(z) in (2) are characterized by the presence of poles at the roots of the equation 

■^w- wgy-fl (6) 

me (o 2 [U 2 (z)-yl\ v ; 

where U{z)= 1 and Jl — ^h ' k/o>. 

OJ 

Hence the solutions of (2) have regular singularities at the roots of (6), which for v = occur within 
the layer unless w << (Oh and the layer is very sharp. Behavior of the solutions of (2) depends 
on the number, location, and indices of these singularities. Any method developed for the study 
of the problem must preserve these basic features of the solutions of (2). 

All the physical properties of inhomogeneous magneto-ionic media are cyclotonic dyads. 
Hence the matrices representing them cannot be diagonalized with all real elements. In other 
words, linear combinations of the type [E x + ip(z)E y ] or [a\E x -\- b\ H x + aiEy + bzHy\ obeying sep- 
arate (decoupled) wave equations do not exist for inhomogeneous magneto-ionic media. 
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Forsterling [1942], Rydbeck [1950], and lately Clemmow and Heading [1954], tried to deter- 
mine suitable linear combinations of the horizontal components of the field vectors for which the 
coupling function would remain small throughout the layer to justify the application of perturbation 
methods. It can be shown that the coupling functions of these equations satisfied by the specific 
linear combinations suggested by these workers have precisely the same number and types of 
poles as L(z) in (2) has. If it were not so, the solutions of the derived equations would differ radi- 
cally from those of the original Maxwell-Lorentz equations. 

A second-order linear differential equation is basically equivalent to a pair of coupled first- 
order equations; moreover, irrespective of the nature of the coefficient functions, a first-order 
equation can always be expressed in closed integral form, while a second-order equation cannot 
be so easily integrated unless the coefficients are elementary functions of the independent variable. 
Hence the methods adopted by these workers amount, at least in principle to, the evaluation of 
the solutions, in closed forms, for any linear differential equation. 

Besides the above-mentioned attempts and the recently developed methods of numerical 
integration [Budden, 1955; Barron and Budden, 1959], this complicated problem has not, as yet, 
been tackled by any worker. In part II of this paper the pair of equations (2) for a parabolic 
electron layer has been investigated, and theoretical expressions for the four observables mentioned 
earlier have been evaluated from the relevant asymptotic solutions. 

2. Wave Equations for the Transverse Modes 

Consider the propagation of a continuous train of em waves of frequency/— co/277 in a magneto- 
ionic medium with the following assumptions: 

(i) The propagation is along the vertical z-axis; 
(ii) The steady magnetic field H P is along the x-axis; 
(iii) Noncoherent scattering is absent; 
(iv) Collision damping v is independent of z; and 
(v) The electron-density N(z) is of the form 



N(z) = N m \l -j\ for - / ^ z ^ /, (7) 

= 0forz 2 >/ 2 . 



It can be shown that (1), under the above set of conditions, reduces to 

a 2 + b 



dz 2 



2 



E x = 0, (8 ) 



dz l 



(z 2 - a 2 ) (z 2 - b 2 ) 



(z 2 -z 2 ) 



£y = 0, (9) 



K 2 = 4tt 2 /i 2 // 4 , h\ = l 2 l(l - iZ)\ 2 m , Z = vlw, 

2ttc 
k m = , the wavelength in vacuum corresponding to the maximum plasma-frequency of the layer, 

0) m 



, = e*Nm 
° m me 



, o»^ = [/XoeH e /m] 2 . c = 3.10 8 m/s, 



a 2 = (l-hl + h 2 h 3 )l 2 , b 2 = (l-hl-h 2 h 3 )l 2 , z\ = {\ - h\ + h$P , 



3 (i-«ZK 

1288 



Equation (8) describe the transverse o-mode for parabolic layers. The magnetic field II, 
has no effect on this mode. This equation has been discussed in details by Rydbeck [1942]. 

Equation (9) describing the transverse e-mode represents the simplest type of magneto-ionic 
propagation in a parabolic layer. It is characterized by the pair of regular singularities at z = ±Zo. 
The nature of propagation for — z () ^ z ^ Zq (the magneto-ionic core) is entirely different from that 
in the regions z 2 > z\ (fringes of the layer). Evidently (9) reduces to (8) for cdh — 0. For w 2 ~ u) 2 m 
± (o 2 , solutions of (9) near z — * — / (bottom of the layer) can be represented accurately by the asymp- 
totic forms. 

3. Solutions of the Wave Equation 

A quadratic transformation of the independent variable in (9) results in an equation having 
one elementary singluarity (exponent difference §) at the origin, one regular singularity at z\ and 
an essential one of grade 1 at «>. Equation (9), or the one transformed as indicated above, belongs 
to a class higher than the confluent hypergeometric equations (Whittaker, Weber, and Bessel 
equations) which accommodate one regular and one essential singularity, or the Mathieu class 
having besides the essential one, a pair of elementary singularities. 

Changing the variables from z to £ = (] -z 2 /z 2 } ) and from E,,(z) to F({)= e -n* 2 -^ 2 • E y (z), (9) can 
be written as 



d 2 F 



U-\)i 



^|+[Mfo£-AF=0, 



(10) 



fhere 



io = iKz 2 = 2inh .,( 1 - h\ + Af), 



M=j[a* + b*-z$-\=^[l-hl-h$ 



1 

"4' 



1/2 

\=-(z?-a 2 )(zl-b 2 ) = -ir 2 h 2 l h*(h*-h$. 



Consider the integral 



I 



{vi-tiWijido, 



(11) 



where 0(rj) obeys the equation 



d_ 



t)(1- 1) )|-( 1 )''0)+{-^ 2 + &-?)t)-M}<0 



At} m = O. 



(Ha) 



Solutions of (11a) can be written as 

8(7]) = i)- M ( 1 - ■n)~ M -" 2 e-^oft(Tj), 



0i(i?) = l- 



M(M- 



-go) + A 



= I- 



M-\ 

M{M-\) + h 



■ 17 + higher powers of 17, 



M-\ 



(1 — 7j)+ higher powers of (1 — 17). 
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Since the series 0i(/n) is analytic throughout the finite range of 77, the behavior of #(17) near the 
singularities 0, 1, and °° are as follows: 

6(r)) -» j]- M as 7) -» 0, 

-+(l-. v yM-ii2 as 77^1, 



-^0 



as 17 



Integral (11) satisfies (10), provided the contour C is so chosen that the bilinear concomitant of 
(17 — €) M an d #(17) vanishes at the extremities of C. From the above behavior of #(77) near its singu- 
larities, contours may be selected for different ranges of f such that (11) identically satisfies (10) 
(see table 1). 

TABLE 1. Contours for different ranges of g so that (11) identically 
satisfies (10) 



Range 


Contour C 


Solution 


0<R ( >t^ 1, 

z 2 ^|z3|;_ 

Magneto-ionic 

core 


(i) 17=1 to 77 = 00, R e {o>0 

(ii) 7? = to 77 = — 00, Reljo < 
(iii) 17 = to 17 = £ 
(iv) 17 = £ to i)= 1 


F-M) 
FM) 

FM) 


£>Wli. 

Magneto-ionic 
fringe. 


(v) 17 = 1 to 17 = 00, R e £ > 
(vi) 17 = £ to 77 = — 00, R e g < 

(vii) 17 = \ to 17 = 

(viii) 17 = to 17= 1 


FM) = F,{£) 

FM) 

FM) 



The general solution inside the core can be expressed by selecting the contour C as consisting 
of a line from 77 = to tj = £, a semicircle around 7) — £, and a line from r) = £ to 17 = 1. Carrying 
out the integrations along the different branches of C and introducing the time factor e llot , the gen- 
eral solution inside the magneto-ionic core can be expressed in terms of the original variable as 



E y (z;t)=A (l-£)exp i|<ot- | z 2 } J +P/ 2 (£) exp i{^ + |z 2 } 



(12) 



where / x (l - z 2 /zg) = J a n (l - zW +1 , 



/ 2 ( 2 / 20 )=2Mz/20) 2 ' , + 1 , 



2r 



a = 6o = 1, 



H 



VttHI-M) 



■exp[-i{2itM+K4}], 



oi = (A+|-ft)/2, fc, = [A-(^+l)f ]/(3/2), 



(»+l) (n + 2)*w* - 



» 2 -i fo-o +A + ^-fo 



(n+l)(« + -) 6„ +2 - 



:f,-|) + A+|-f, 



o„ +1 + (M — ra)a„ = 0, 
b„ +l -Zo(M + n+l)bn = 0. 



The right-hand side of (12) represents two wavy disturbances of complex amplitude ratio P, 
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moving in opposite directions with complex phase factors (— Kz 2 /2 — i In /i) and (Kz 2 /2 — i In f 2 ), 
respectively. For very long waves such that u> < con, the location of the singularity at z§ with i> = 
is outside the physical layer of half-width /. Hence for such waves, the composite disturbance 
at the bottom of the layer z = — l can be expressed as 



Ey(-l-t)-- 



i+^^-^PdKP) 



E° y (-l; t), 



(13) 



where E° y ( — l; t) is the state of the wave incident on the layer. If 



Re 



^ _ M-l)exp(iKF 
PfA-D 



where R e and </> e are real quantities, the physical observables are related to the above variables as: 
R e is the amplitude of the downcoming wave at z = — /, relative to that of the upgoing wave at that 
point = R e e x (\ where R e is the reflection coefficient for a collision free layer, and Xe is the total 
absorption in the layer due to collisional damping. If r e be the time delay observed on the CRO 
screen and hf } the height of the bottom layer from the transmitter-receiver, 

2Kb . </>e 

Te = h— , 

C (x) 

where c=3.10 8 m/sec. Hence solution (12) of (9) furnishes information about all the observabies 
of the transverse e-mode, those of the o-mode being derivable from these results of putting ojh = 0. 



4. Asymptotic Solutions 

For waves used in normal radio cummunications, o>// <$ w ~ a) m and hence \zl\// 2 ^ 0. Solu- 
tions of (9) for these frequency bands are given by the asymptotic behavior of E u (z; t). In this 
region of frequencies consider the integral 



/, 



e ^o* v{r)) e/n, 



(14) 



vhere 



va-v)^+[{oV^({o + ^V + (M+2)^+[2{ oV -(^o + A + ^]v = 0. (l 4a ) 



Solutions of (14a) can be expressed as 

v{r)) = 7]-^ M) (l-r?)*'- 1 / 2 e^o ]T y n ^ 

o 

with 70 = 1 , 



7i : 



M 



, and 



(/i + 2)(7i+l-M)y /l+2 + 



(/i+l)(/i + --fo)+A 



7n+i — ngoy n = 0. 



Integral (14) will satisfy (10) if the contour C is properly chosen. Since the real axis in the 
-n-plane is the cut plane of the integral in (14), the solution for the sector ^ arg(— £ £) < tt will be 
different from the solution in the sector it ^ arg (— fof) < 2tt. For (— £of)= |£o£|e Ja and R e M < 0, 
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the contour C may be chosen to consist of 

(a) a line G, joining 77 = to 77 = — \r\e~ ia , 

(b) a line C 2 , joining rj = — |r|e~' a to 17= 1— |r|e~* a , and 

(c) a loop C;i, in the positive direction round r\ = 1 starting and ending at 77 = 1 — |y|e~ ia . 

As |y| — > °°, the contribution from Cz will vanish, and the solution of (10) for f — > — 00 can be ex- 
pressed as 



F(f) 



-/; 



e -i?l ^ . ~-(l+M) 



£ (1 - ijf-^erto ' 7nT?" 



6/17 



+ 



f(l + ) 
•^ 1-ooe" 



e -r,^ 1 _ T? )A/-l/2 



2 77- (1+M) e^oy n T7 w 



G?T7* 



For |f I > 1, the asymptotic solution can be obtained by putting 77 = (ai/fof) in the first integral and 
l—r) = (—U2l£o£) in the second. If further, (1/fof) is written as e i7r /(— f f) for < a ^ 77 and as 
e~ ln l(— f f ) for 7r < o: ^ 27T, the form of the solution as z —> — I can be expressed as 



em ~ n- M) • (- totr ■ Vi(- f of) • ex P • [y 2 f of ] 



where 



+ 



2tt 



r(-M+y 2 ) 



jtt • (- fof)~ M ^ /2 r 2 (- f of) exp • [- Yaf of ] , 



V( ,e. . , A M(M-1) . ■ ■ (M-b+1) 
F l( -6f)=l + S ^^ «». 

V 2 (~ fof ) - 1 + 2 (" D" ^jj^y, ft., 



and 



2a n x n s £ y n [^n(l_^)M-l/2 . e xf ] ? 



£/3n* n = £ Ml -^) n " 1 - M ' e* 1 "^]. 



Introducing the time factor e i(ot and passing on to the original variables, the form of the dis- 
turbance as z— > — / for q)h ^ co ~ co m can be expressed as 



where 



E y (z; t)-f 3 (z) exp - 



MO)£- 



£z 2 



+/ 4 (z)exp. i\ 



i\ a)t + 



Kz 2 



y(-m + - 



f 4 (z) ~ Y(-M) [iK(z 2 -zm M exp 



.W +irM _ \-iMW-€o) 

i 2 ^ 7rm K(z*-z%) 



m^ a K-mjM-m 



■nM 



K(z*-zl) 



iM=M +1/4. 

Since 2 < near the bottom of the layer, the first term represents the upgoing (incident) wave 
while the other represents the downcoming or reflected wave. Hence if E\(— I; t) represents the 
state of the incident wave at z——l, the composite disturbance at that point can be written as 



E y (-h t)= [l+Ree-^e] £•«(-/; t), 
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(15) 



yher 



R e e-*e=\f 4 (z)exp(iKz 2 )lf 3 (z)\ i 



(15a) 



For the o-mode, given by (8), similar expressions can be deduced from above by letting co//— » 0. 
If R$ and R 2 be the reflection coefficients of the two modes in the absence of collisional damp- 
ing, it can be shown that 



77/co 



In 



In 


r m 




R 2 e 1 


d 


I- 


~R 2 e\ 





=/o=l- 



af 



ix> 



D2 ~f e ~^ . .2 






(16a) 
(16b) 



"m w m 



Here R 2 and /?^ are unity throughout the whole range of frequencies except at the critically penetrating 
bands (co 2 ~ co 2 for o-mode and co 2 ~ o£ — w 2 H for the e-mode), where they suddenly fall to half the 

value. Variations of the left-hand expressions in (16a) and (16b) with hi=-^- have been shown in 
figure 1. 

The total absorptions x and x P for the two modes are given by 

X = ln (RoIRo) and x e = ln (£/Jt.). 

Separating the right-hand side of (15a) into real and imaginary parts and retaining only the dominant 
terms, it can be shown that 



Co 2 


[{- 


«,-? 


[{' 


8= ^ 


I 

c 



l-fln(i)- 

^ J vo/ 



2 



(17a) 
(17b) 



Figure 2 shows the variations of Xo and \e with hi = a)/(o m (for v = 0). 

If To and T e be the observed time-delays between the first reflected waves of the two modes 
and the ground pulse it can be shown that 



To 



2h h 

— +g 


M- 


-!»(£)] 


2hb . 


M- 


-I-©] 



77 
CO 

77 
CO 



(18a) 



(18b) 



where hb = height of the bottom of the layer. Variations of the function y= \t h — with 

c 2a> 

with h 2 = — (for ^ = 0) for the two modes have been shown in figure 3. 
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FIGURE 1. Reflection coefficient as a function of 
frequency. 



FIGURE 2. Absorption coefficient as a function of 
frequency. 
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FIGURE 3. Time-delay as a function of frequency. 
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Expressions (16a) to (18b) furnish informations about all the quantities observed in the radio 
exploration of the ionosphere. Though the above analysis has been carried out only for a para- 
bolic layer, the method can be extended to any other type of electron layer. 

5. Appendix. List of Symbols and Abbreviations 

Electric vector of the wave E (x, y, z; t). 

Permittivity of free space eo- 

Permeability of free space jxq. 

Steady magnetic field H e . 

Collision frequency v. 

Electronic mass and charge m y e. 

Frequency of the traveling wave. /= co/27T. 

Velocity of light c. 

Electron density N. 

Maximum electron density N m . 

Half- width of the layer /. 

Height of the bottom of the layer A&. 

Reflection coefficients 7? 2 , R 2 . 

Total absorption X > X e - 

Observed time delay on the CRO screen to, r e . 

Ratio of the axes of the polarization ellipse $ft . 

Tilt of the polarization ellipse ^p. 

Abbreviations: 

(t>l l = N m e 2 lm€ , \ m = 2ircla) m , oj h = |e/icH e /ra|, 

z=\- ivfa , k* = why i 4 = fer) 2 /(i - iz\ 

h\ = / 2 /(l - iZ)\l, hi = (1 - iZ)a> 2 /(*) 2 m, hi = oA/(l - iZ)a) 2 m, 
a 2 = (l-ht + h 2 h 3 )l 2 , b 2 = (l-h 2 - h 2 h 3 )l 2 , z% = (1 - h\ + A 2 )/ 2 , 
M= l ^(l-h 2 2 -hl)-\ = iM-\, 

A = -7T 2 h 2 hl(h 2 2 -hl), 

£o=iKzl, g=llch 2 (for *> = 0), 

r , CO 2 (O 2 CO 2 
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